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In 2011 and 2012, three reactor neutrino experiments. Double Chooz, Daya Bay and RENO 
showed positive signals of reactor neutrino disappearance and measured a mixing parameter sin^ 2^13 
at average baselines 1.05, 1.65 and 1.44 km, respectively. It is possible to measure effective Amli 
(Am^ defined in two flavor oscillation formula, hereafter referenced as Am|i) from distortion of 
neutrino energy spectrum {E dependence of the oscillation) in those experiments. However, since 
it requires a precise energy calibration, such measurements have not been reported yet. Am|i can 
also be measured from baseline (L) dependence of the neutrino oscillation. In this paper Am|i 
is measured from disappearance probabilities of the three reactor experiments which have different 
baselines, to be x 10"^ eV^^ where the errors are two (one) dimensional uncertainties. 

This is consistent with Am|2 measured by z/^ disappearance in accelerator experiments. Importance 
of Am|i measurement and future possibilities are also discussed. 



I. INTRODUCTION 

Neutrino oscillation is, so far, the only firm phe- 
nomenon which is not accounted for by the standard 
model of elementary particles, which assume neutrinos 
as massless. The neutrino oscillation is, like other os- 
cillations, such as O ^0 , 5° ^ ^(CP violation), 
d ^ s (Cabbibo angle), {ud) <^ {du) (isospin), B <^ W^, 
(Weinberg angle), p{^)e{\) <^p(i)e(t) in hydrogen atom 
(21 cm HI line) etc, assumed also to carry very important 
physics and we should be able to learn much about our 
world from it. 

There are six parameters in standard three flavor neu- 
trino oscillation [1]. Three mixing angles between flavor 
eigenstates and mass eigenstates: 6^12, ^13, 6^23, one CP 
violating imaginary phase: 5 and 2 independent squared 
mass differences: ^m'^i. = m'j — m^, where are neu- 
trino masses (mi, m2, ms) of the three mass eignes- 
tates (i^i, ^3) which correspond to the largest compo- 
nent of {ye^h>^^h>^)^ respectively. Before 2011, sin^26>i2, 
sin^ 2^23 5 ^^21 5 I ^^32 1 were measured by various ex- 
periments and observations [1]. 6^13 was known to be 
small, sin^ 2^13 < 0.15, from Chooz reactor neutrino ex- 
periment [2]. In order to measure J, to determine mass 
hierarchy and to solve ^23 degeneracy, 6^13 has to be rel- 
atively large. Thus finite value of ^13 had been eagerly 
sought for. 

The years 2011 and 2012 will be regarded as an epoch 
making ones of neutrino experiments. T2K group showed 
^ ^ appearance candidates over 1.5 expected 
backgrounds in June 2011 [3]. MINOS group showed 
also an indication of appearance [4]. Double 

Chooz showed an indication of the reactor neutrino dis- 
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appearance in November, 2011 [5j. Daya Bay and RENO 
showed more precise disappearances on March and April 
2012, respectively 0,0. In June 2012, at neutrino con- 
ference held in Kyoto, Double Chooz [8], Daya Bay [9[, 
T2K [10] and MINOS [IT] updated their sin^ 2l9i3 mea- 
surements. All those results show relatively large Ois and 
have opened up a brilliant path to future neutrino exper- 
iments. 

Reactor neutrino oscillation probability is expressed as 
follow. 

PniPe ^ Pe) ~ 1 - sin^ 2^13 sin^ ^^L, (1) 

where, Ej^ is neutrino energy (~ a few MeV) and L is 
baseline (1 ~ 2 km). Amg;^ "effective" squared 

mass difference often mentioned as Am^^^ in two flavor 
analyses. In three flavor oscillations, it is an average of 
Am^i and will be described in the next section. 

In principle both sin^ 26>i3 and Amg^, can be measured 
from the oscillation. However, all the reactor neutrino 
experiments use Am32, which was measured by dis- 
appearance by MINOS group [12], as Am^^ to extract 
sin^ 2^13 since the difference between them is in an or- 
der of |Am2i/Am32| ~ 3% which is smaller than current 
accuracy of the measurements. 

It is important to measure Am^^ independently from 
Am32 to check if the standard three flavor scheme is cor- 
rect. If there is significant difference between Am32 and 
Amli, it indicates an existence of new physics. More- 
over if they are measured with precision of 1% or better, 
mass hierarchy can be solved and cos S may be measured 
as described in the next section. 

Experimentally, Arfi'^i can be measured by analyzing 
either E dependence of the oscillation or L dependence 
of the oscillation. Both methods use independent infor- 
mation, namely energy distortion and normalization, so 
that combining both analyses, accuracy of the Amg^ will 
improve. The former requires a precise energy calibra- 
tion and no results have been reported yet. The latter 
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analysis can be performed by combining currently avail- 
able disappearance information at different baselines. 

In this paper Arhl^ is measured using reported 
sin^ 2^13 and baseline of each reactor experiment. The 
contents of this paper is based on our poster presenta- 
tion in the conference of neutrino 2012 [13|. In next sec- 
tion, neutrino oscillation formula is described stressing 
on relation between Am'^i and importance 
of Arhli measurement will be discussed. In section- 
Ill, how we treat reactor neutrino will be described. In 
section-IV, most recent Double Chooz, Daya Bay and 



RENO results 0-[9j are combined and Am^^ is extracted. 
In section-V, a new experiment is proposed to measure 
Am^^ more precisely by using the reactor complementar- 
ity. 

II. NEUTRINO OSCILLATION FORMULA 
AND EFFECTIVE Am^ 

The mixing matrix between flavor eigenstats and mass 
eigenstats is [l,]. 



C12C13 S12C13 5136 



Uaj = -S12C23 - Ci2S23Si3e'^ C12C23 - Si2523Si3e'^ S23C13 , (2) 



512523 - Ci2C23Si3e*^ -^23^12 " 5i2C235l3e*^ C23C13 



where a is an index of flavor (a = e, /i, r) and j is an 
index of mass eigenstates (j = 1,2,3), cjk = cosOjk^ 
Sjk = sin Ojk^ and tjk = tdiii Ojk will be used later. S 
is so called CP violating imaginary phase. Currently 
these parameters are measured as follows O12 ~ 

34°, 6>23 ^ 39°, 6>i3 - 9°, Am^i - 7.5 x IQ-^eV^ and 
lAm^- 2.4 X lO-^eV^ 

Neutrino oscillation probability going to the same fla- 
vor is expressed by following formula. 



P{Ua ^ i^a) = 1 



j>k 



A,, 



jk 



(3) 



where Ajk = The second term in right hand side 

is called disappearance probability. This oscillation for- 
mula is valid for both neutrino and antineutrino cases. 

Reactor neutrino experiments use Ve generated by j3- 
decays of the fission elements in the reactor core. Energy 
of the neutrino is a few MeV. At around the first os- 
cillation maximum of A32, survival probability of Ve is 
expressed as. 



1 - sin2 2^13 [cl, W ^ + si, W ^ j + 0(10-3). 

(4) 

On the other hand, the survival probability of high energy 
which is produced by accelerator is. 



(si2 + ^13^23 sin 26>i2 COS 5) sin^ ^ 
+(c^2 - 513^23 sin 26>i2 cos S) sin^ 



0(10- 



(5) 



Usually oscillation data are analyzed by assuming two 
flavor oscillation formula. 



P{j^a ^ ^a) = 1 - sin^ 26>sin' 



Am^L 

4.E ' 



(6) 



and the measured mass square difference corresponds to 
a weighted mean of |Am32| and lAmg;^! [l^]^ 



Am3i |Am3i| 



^2 



I Am 



32 I 



Am^2 



{si. 



-§13^23 sin 26^12 cos S) | Amli | 
+ (ci2 — 513^23 sin 2^12 cos S) I Am32 1 



(7) 



They are called effective Aw?. Note that Afh? is not 
a difference of the mass square and is positive definite. 
Since there is a relation 



Amii 



Am32 



Am|i, 



(8) 



in the standard three flavor scheme, the difference of 
Amg^L Am32 is expressed as follows. 



2(Amii 



Ami2) 



Affile 



Arfil^ 



2 cos 2i9i2|Am|i| 
|Amii| + |Ami2| 



± (1 — 513^23 tan 2^12 cos 5) 



±0.012 X (1±0.3), 



(9) 



where the overall sign depends on mass hierarchy. If 
Am^i > Am32, it is normal hierarchy, and vise versa. 
In order to distinguish the mass hierarchy cases, it is 
necessary to distinguish the separation of 1.7^3.1% de- 
pending on S. Am'^2 been measured with precision of 
~ 3.5% [14]. So far there has been no reported measure- 
ment of Am'li and this paper is on the first measurement 
of it. If difference between Am'^i and Am32 is larger 
than 1.6%, it can not be explained by the standard three 
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"^31 



and Am32 are 



flavor oscillation scheme. If both A 
measured with accuracy 1% or better in the future, the 
mass hierarchy and cos S can be measured. 



III. REACTOR NEUTRINO OSCILLATION 

In nuclear reactors uranium and plutonium perform 
fission reaction; after absorbing a thermal neutron they 
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break up into two large nuclei called fission products, and 
two or three neutrons which sustains the chain reaction 
of the fission. The fission products are generally neutron 
rich nuclei and unstable. They perform ~ 6 /3-decays on 
average before becoming stable. In each /3-decay, a 
is produced. On the other hand, 200 MeV of energy is 
released per fission [16], which means ~ 6 x lO^^i^e are 
produced every second in a typical power reactor with 
3 GW thermal energy. At 1 km from such reactors, v^. 
fiux amounts to ^ 5 x 10^/s/cm^. The energy spectrum 
of the reactor neutrinos is a sum of the energy spectrum 
of neutrinos originated from the four fissile elements. 



E 



(10) 



-235U 238U 239py 241py 



where fi[Ey) is reactor neutrino spectrum per fission 
from fissile element i and is a fraction of fission 
rate of fissile element i. There is a relation = 
1. For equilibrium light water reactors, /3i are simi- 
lar and we use the values in Bugey paper [16|, namely 

235U : 238U : 239p^ . 241p^ ^ q 533. q Q^g. q 323. 

0.056. In this study, fi{E) is approximated as an expo- 
nential of a polynomial function which is defined in , 



(11) 



In reactor neutrino experiments, usually organic liquid 
scintillator is used to detect Uq. It is rich in free protons 
and reactor fj^ performs inverse /3-decay interaction with 
a proton. 



- p ^ 



(12) 



This is an inverse process of neutron /3-decay (IBD) with 
very small q'^ and the cross section is precisely calculated 
from the neutron lifetime [18]. In this analysis, infor- 
mation of absolute normalization is not necessary. The 
energy dependence of the IBD cross section is. 



(Tibb{E^) cc (^,[MeV] - 1.29)y^ - ^ lA. 

(13) 

The disappearance probability,^, can be related to the 
oscillation parameters such as |19| . 



Pd = sin^ 20 



J S,{E)aiBB{E) sin^ (^) dE 



J S,{E)aiBB{E)dE 
= sin2 2i9A(Am^L). (14) 

Reactor measurement of Ois is a pure sin^ 2^13 mea- 
surement in contrast to accelerator based measurements 
which depend on unknown parameters. It means that by 
combining reactor results and accelerator results, infor- 
mation of such unknown parameters can be derived [20] . 
Under such motivations, several reactor-^13 experiments 
were proposed in the past [2l| and now Double Chooz, 



Daya Bay and RENO experiments have published pos- 
itive results of the reactor neutrino disappearance and 
measured sin^ 26^13. These experiments make use of the 
same concept to reduce systematic uncertainties signifi- 
cantly [22] over the previous experiments of Chooz [2] and 
Palo Verde j23|. That is, they construct far detector(s) 
at around oscillation maximum and measure the "oscil- 
lated" spectrum. On the other hand, near detector(s) 
with same structure as the far detector is constructed at 
around a few hundreds of meters from their reactors to 
measure the neutrino spectrum before the oscillation. By 
comparing the data taken by the near and far detectors, 
the only effect caused by the oscillation can be derived 
by canceling systematic uncertainties of reactor neutrino 
fiux and detection efficiencies. The ffux- weighted average 
baselines of the far detectors, (L), for the three reactor 
experiments are 1.05 km for Double Chooz, 1.44 km for 
RENO and 1.65 km for Daya Bay, respectively. The av- 
erage baseline of RENO experiment was calculated using 
neutrino ffux shown in and distances between the far 
detector and each reactor. Published values are used for 
Daya Bay and Double Chooz. 



IV. COMBINATION OF THE RESULTS FROM 
THE THREE REACTOR EXPERIMENTS 

From measured disappearance probability, an allowed 
line can be drawn in sin^ 20 — Am^ parameter space us- 
ing the relation (p!4|) . Since the baselines are different for 
the three reactor neutrino experiments, there are three 
different allowed lines as shown in fig.-ITJa). The point of 
intersection indicates the solution of Am^^ and sin^ 2^13. 
In real experiments, due to errors, the three lines do not 
cross at same point. When combining different reactor 
results, values are calculated by using following for- 
mula for each point of the parameter space. 



X 



E 

fe=exp. 



sin^ 2eA{Am^{L)k)-{Pd)k 



(15) 



where k is index of the three experiments and dk is mea- 
surement error of experiment- /c. Fig.-[Hb) shows contour 
of the significance in case each experiment measures the 
disappearance with 0.5% accuracy. In this case, Am^^ 
can be determined with precision of ~ 23(9)% with two 
(one) dimensional uncertainty. The large difference be- 
tween the one dimensional error and the two dimensional 
error is because the shape of one a contour island has long 
tail as shown in fig.-[Hb). Since there are two parameters 
to measure, at least three experiments are necessary to 
redundantly measure the parameters. 

In the actual analysis, Pd is not directly written in pa- 
pers and it is calculated from measured sin^ 2^13 and fiux- 
weighted mean distance (L). In their papers, sin^ 2^13 
were derived by assuming the MINOS Am32 [12|. Re- 
lations between these parameters and the disappearance 
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FIG. 1. Hypothetical sensitivity plots assuming true parame- 
ter values are Am^ = 2.32 x 10~^ eV^ and sin^ 26 = 0.092. (a) 
Allowed lines if disappearance probabilities are measured as 
expected, (b) Allowed regions after combining the three reac- 
tor experiments. The contour lines correspond to, from inner 
to outer, Icr, 2a and 3a significances. Disappearance error of 
0.5 % is assumed for each experiment. Amii is expected to 
be measured with ~ 23(9) % accuracy corresponding to two 
(one) dimensional uncertainty. 



probability, and allowed line are shown below. 

Pd = sin2 2^13A(Am^2WFar) 

= sin^ 2i9A(Am2(L)Far). 



(16) 



The calculated disappearance probabilities are shown in 
table-U together with other parameters. 



TABLE I. Parameters of the three reactor neutrino experi- 
ments. 



Item 


Double Chooz 


Daya Bay 


RENO 


(L) Far [km] 

sin^ 20i3 

Pd@(L)Far 


1.05 
0.109 ±0.039 
5.5 ±2.0% 


1.65 
0.089 ±0.011 
7.0 ± 0.9% 


1.44 
0.113 ±0.023 
8.2 ± 1.6% 



FIG. 2. Same as fig.-[T]but disappearance probabilities and 
their errors are calculated from measured sin^ 2^13. 



rameters and their errors are, 



Ami, = 2.99i}i3( 



+0.86 
0.88 



sm' 



X 10" 
0.013)' 



+0.014^ 



(17) 



where the errors are for two (one) dimensional uncer- 
tainty. This result is consistent with within one a 
and there is no deviation from 3 flavor oscillation within 
this accuracy. The most probable sin^ 2^13 value coin- 
cides with the Daya Bay result but this sin^ 2^13 has 
meaning that it was derived without assuming Am3;L- 
The minimum is 0.43 with one degree of freedom 
which means that the results of three reactor experiments 
are consistent with each other. 

Baseline dependence of observed disappearance prob- 
ability and various expectation lines are shown in fig.-lSl 
This figure clearly shows the relation of the calculated 
disappearance probabilities and expected oscillation pat- 
terns. The meaning of disappearance probability is also 
described in its caption. In the near future, errors of the 
experiments are expected to improve much and the oscil- 
lation pattern will be determined much more precisely. 



The sin^ 2^13 were measured using both near and far 
detector at Daya Bay and RENO but only far detector 
was used in Double Chooz experiment. It is important to 
point out that although the reactor experiments assume 
Am32 to extract sin^26>i3, this analysis is independent 
of the assumption for the first order. If the experiments 
used different Am^ , they would obtain different sin^ 2^13 
but the Pd calculated by the equation (p!6|) would be the 
same. Am32 was used just as a reference point. Fig.-|2] 
shows the combination of the three reactor experiments 
calculated this way. The most probable oscillation pa- 



V. FUTURE POSSIBILITIES 

It is important to evaluate how precisely we can mea- 
sure Arfi'^i since it may resolve the mass hierarchy com- 
paring with Am32 in the future. In order to make the 
most of the reactor complementarity, we studied a case 
to add a fourth experiment and calculated an optimum 
baseline to measure Amg^ by combining with the cur- 
rent three experiments. Fig. -H] shows dependence of the 
two dimensional uncertainty on baseline of the fourth 
experiment. The accuracy improves rapidly when L ex- 
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FIG. 3. Baseline dependence of Ve survival probabilities. 
Dashed and/or dotted lines are expected oscillation pattern 
calculated using sin^ 2^13 measured in each experiment and 
MINOS Ami2- The calculated disappearance probabilities 
correspond to the values of the expected lines at their flux- 
averaged baselines. The solid line is expectation from the 
most probable Am|i and sin^ 2^13 measured by this analysis. 
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FIG. 4. Dependence of the two dimensional uncertainty on 
baseline of fourth experiment. 



ceeds Daya Bay baselines of 1.6km and reaches to '^7% 
at 2.5 km. This is because that the tail of the island 
in the sensitivity contour plot vanishes thanks to the 
almost perpendicular intersection of the fourth allowed 
line. Fig.-|5] shows sensitivities with the fourth experi- 
ment with baseline 2.5 km. Since the combined allowed 
region no more has tails, one and two dimensional errors 
become similar. This baseline is mere 1.5 times of Daya 
Bay far detector baselines and thus it is not unrealistic. 
After detailed energy calibrations are finished, the cur- 
rent reactor experiments will perform spectrum shape 
analysis to extract Am^]^. Since the shape analysis and 
the baseline analysis use independent information, the 
accuracy of the Amg^ is expected to improve by com- 
bining them. However, in order to determine the mass 
hierarchy, one step more improvement of the accuracy 
will be needed for both Amoi and Amoo. 



FIG. 5. Same as fig.-l but a hypothetical fourth experiment 
with baseline 2.5 km is added. Amli can be measured with 
precision ^7%. 



VI. SUMMARY 

In this paper, measurements of Am^^ by using the 
baseline differences between currently running reactor 
neutrino experiments were studied and we obtained fol- 
lowing results. 

(1) About motivations, independent measurement of 
Am3;L is important with following reasons, (i) The stan- 
dard three flavor oscillation scheme can be tested, (ii) 
Consistency among results from reactor neutrino exper- 
iments can be checked. In order to perform the con- 
sistency check, at least three experiments are necessary, 
(iii) It may resolve mass hierarchy and give information 
of cos(^ in future experiments, (iv) The reactor com- 
plementarity method uses independent information from 
spectrum shape analysis and accuracy of Amg^ will im- 
prove if they are combined. 

(2) The current data from Daya Bay, RENO 
and Double Chooz were combined and Am^^ = 
2-99tl;58(-o;88) was obtained, where the first (sec- 
ond) error was two (one) dimensional uncertainty. This 
is consistent with Am32 measured by the accelerator ex- 
periment. Minimum was small which shows that the 
results of the three experiments were consistent with each 
other. 

(3) It was demonstrated that Am^^^ could be measured 
with precision ~23(9) % by combining the three exper- 
iments if the error for their disappearance probabilities 
would reach 0.5%. 

(4) A new detector with baseline 2.5 km will reduce 
both the one and two dimensional uncertainties of Am^]^ 
to ^ 7% by combining with the current three reactor 
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experiments. 
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